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Abstract. On the basis of a hydrodynamical model analogous to that in critical 
fluids, we investigate the influences of shear flow upon the electrostatic contribution 
to the viscosity of binary electrolyte solutions in the Debye-Huckel approximation. 
Within the linear-response theory, we reproduce the classical limiting law that the 
excess viscosity is proportional to the square root of the concentration of the electrolyte. 
C$ , We also extend this result for finite shear. An analytic expression of the anisotropic 

structure factor of the charge density under shear is obtained, and its deformation 
' ■ at large shear rates is discussed. A non-Newtonian effect caused by deformations of 

the ionic atmosphere is also elucidated for 7737 > 1. This finding concludes that the 
Q ' maximum shear stress that the ionic atmosphere can support is proportional to A^ 3 , 

where 7, A^ and 775 — \ 2 D /D are, respectively, the shear rate, the Debye screening 
length and the Debye relaxation time with D being the relative diffusivity at the infinite 
dilution limit of the electrolyte. 
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1. Introduction 

It is deeply recognized that the long-ranged Coulomb interaction between ions gives 
rise to the viscosity enhancement of electrolyte solutions [H 12 EH E] ■ When a velocity 
gradient in the solution will deform an ionic atmosphere (a cloud of couterions around 
an ion), which would otherwise possesses a spherical symmetry, electrostatic forces as 
well as thermal motion tend to restore the atmosphere to its original form. Because of 
a finite relaxation time of the ionic atmosphere, these two effects can be balanced in a 
stationary state. This effect is closely related to the relaxation effect in conductance, and 
also it is sometimes called the "electroviscous" effect by analogy with the electroviscous 
effects found in charged colloidal suspensions In the case of very low concentration, 
such an ionic effect is proportional to the square root of the concentration, as was found 
many years ago by Falkenhagen |Hj and Onsager and Fuoss |Sl Ej theoretically, and by 
Jones and Dole experimentally jHj. 

The concentration dependence of the transport properties (such as the conductivity 
and the viscosity) of electrolyte and polyelectrolyte solutions still offer significant 
challenges. Today, there exists a large number of theoretical and numerical studies 
that take into account realistic structural effects of ions and polyions. Among them, 
one of the most celebrated is the mode coupling theory combined with accurate pair 
correlation functions of liquids .9] , which has recently made considerable progress in this 
area [TQ1 HH 112 EH Ej- With regard to the study of viscosity, however, most of the 
literature has concentrated on the intrinsic viscosity. 

To the best of our knowledge, nonlinear responses of electrolyte solutions in velocity 
field remains largely unexplored so far, even in the infinite dilution limit. The aim of this 
paper is therefore to study how a strong shear flow modifies the "electroviscous effect" 
mentioned above. Because analytic calculations of general cases are impossible, we will 
consider in this paper a binary electrolyte solution in very low concentration, for which 
analytic expressions of various statistical quantities are available. Our formulation is 
in principle applicable over an unrestricted range of field intensities as far as the local- 
equilibrium condition is not severely violated. Although the present work is limited to 
the specific simple system and thus might be less noticeable in the practical point of view, 
we expect that this work may prompt further studies towards better understandings 
of rheological properties in more important systems such as concentrated electrolyte 
solutions, polyelectrolyte solutions [THj and charged colloidal suspensions [TBI E! under 
shear. 

This paper is organized as follows. In Sec. II, we set up Langevin equations 
for fluctuating concentration of ions combined with the free energy functional which 
accounts for electrostatic interactions within Debye-Hiickel approximation. Our model 
is directly employed from that in critical fluid binary mixtures |18[ ITU] . The linear 
response is first studied in Sec. III. We show a smart derivation of Falkenhagen- Onsager- 
Fuoss limiting law jHl El El for binary electrolyte solutions. In Sec IV, an anisotropic 
structure factor of charge density and the excess viscosity are explicitly obtained as a 
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function of the shear rate. The limitation and the experimental relevancy of our results 
are discussed in Sec V, followed by our brief summary in the last section. 

2. Model Equations 

2.1. Free Energy 

Consider an electrolyte solution consisting of only two kinds of ions dissolved in a 
continuum solvent of dielectric constant e, each of which carries charges Ze and — e 
(e is the elementally charge). The average number density of the two ion species are 
written as (n + ) = n and (n_) = Zn, respectively. Then the total ion number density 
is given by n to t = (Z + l)n. The free energy of this system is composed of the entropic 
and electrostatic contributions: 

f + lj Pc*] dr. (1) 

The first term of the right-hand side of Eq. is the entropic contribution of ideal 
mixing given by jTH] 

/ = k B T^ [n a \n(n a a 3 ) - n a , (2) 

a 

where n a (r) is the number density of the ion species a, T the temperature, k B the 
Boltzmann constant, and the two kinds of ions have a common radius a. The charge 
density is thus given by p c ( r ) — Zen + (r) — en_(r) = eZnif>(r), where we have introduced 
the variable if) as the order parameter of the charge density fluctuation. The second term 
in the free energy F is the electrostatic contribution. The electrostatic potential satisfies 
the Poisson equation 

Ait 

V 2 0(r,t) = -— Pe (r,t). (3) 

Note that the time dependence of the electrostatic potential entirely comes from slow 
variations of the charge density p c (r, t) due to relative diffusion between oppositely 
charged molecules. Introducing the Coulomb operator v(r) which satisfies V 2 t>(r) = 
— 47r<5(r), we can write the formal solution of Eq. (JHJ) as 

0(r, t) = -J dr'v(r - r')p c (r / , t). (4) 

Assuming the total number density of ions is nearly constant everywhere: n + (r) + 
n_(r) = ntot, Eq. (J2J) can be expanded in powers of if). Up to the quadratic order, we 
find / = ip 2 /2x, where \ = (Z + is the osmotic compressibility. Substituting 

these equations into Eq. ((TJ), we obtain the free energy relevant to our purpose in the 
form 

(3F = J dr^ip 2 + Z2 ^ n2 J dv J dr>(i>(r - r')^(r'), (5) 

where /3 = l/k B T and the Bjerrum length £ B = e 2 /eksT is defined as the length at 
which the Coulomb energy between two elementary charges becomes comparable to 
thermal energy k B T . 
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2.2. Dynamic Equations 

To describe the dynamics of the charge fluctuation 5p c = Zen5ip and to calculate the 
shear viscosity under steady shear flow, we start with the following generalized Langevin 
equations (see also Appendix) 

= _ v . V ^-LV 2 ^(/3F)+£, (6) 

P-Qf = "Vp - k B Tif;V—(f3F) + 77o V 2 v + f , (7) 

where the generalized diffusion equation for if> is now coupled with the velocity field of 
the fluid v through the reversible term — v- Vip. The pressure p in Eq. (JZj) is determined 
so as to satisfy the incompressibility condition V ■ v = 0. The Gauss-Markov thermal 
noise source 9(r, t) and f (r, t) are related to the kinetic coefficients L and the zero-shear 
viscosity r] via the usual fluctuation-dissipation relations given by 

(e(r,t)e(r',t')) =-2LV 2 5(r-r')5(t-t'), (8) 
(/ i (r,t)/ i (r',0) = -2vok B T\/ 2 5(r - r')6(t - t')5 iy (9) 

The second term on the right-hand side of Eq. (jJJ) represents the force exerted to the 
fluid by the deviation of if) from its equilibrium value (if>) = 0. This term should be 
incorporated in Eq. (JJJ) as the counterpart of the reversible mode-coupling term — v • Vif) 
in Eq. (jHJ) in order to ensure that the equilibrium distribution of if) and v is given by the 
Boltzmann distribution exp(— F//cgT — §/ <irv 2 //cBT). Neglecting the random force, 
we can rewrite Eq. (J7|) in the form which suggests the momentum conservation 

J^(pv) + v-n = o, (io) 

where the stress tensor is written as 

k B T 2 \ e 

The second term —(e/4n)di<f)dj<j) comes from the long-ranged interionic Coulomb 
interactions and is nothing but the Maxwell stress tensor [20J. When the average velocity 
field (v) = 7ye x is imposed in this system, the corresponding shear stress is thus given 
by 

/tt \ t I dtp d(j)\ 

°* = -<M = w + to\tedjj/- (12) 

Note that the second term in this equation gives rise to the electroviscous effect in a 
sheared electrolyte solution. 



U ij = \ P+ -^—^ 2 S ij ~ ~ VoidiVj + djVi). (11) 



3. Linear Response Theory: The Limiting Law 



In this section, we shall discuss the electrostatic contribution to the viscosity within the 
linear response theory. Because we are interested in the macroscopically homogeneous 
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state ((■0) = 0) far from its critical point, we can linearize Eq. to obtain 

^ = -u • V0 - D(V 2 + 4)0 + 9, (13) 

where u = (v), and where we have defined the relative diffusivity D = L/x and the 
inverse Debye screening length kd = 1/Xd as 

2 A o 72-2 4nZe 2 n tot 

k d = Ani B Z n X = j — 7f — • (14) 

ek B T 

In the absence of the macroscopic flow u = 0, Eq. (fT^j) can be easily solved to yield the 
dynamic structure factor 

Here, the Fourier transform of ip is defined by 

V>M) = J ^exp(zk-r)^ k (t). (16) 

The intrinsic viscosity including the electrostatic contribution 77 can be calculated by 
making use of the Green-Kubo formula [THj : 



1 



00 



71 = r]o+ k B TJo dt J dr \ n * y(r ' t)n ^ (0 ' 0) ) ' (17) 

where the electrostatic shear stress is given by fl xy = e/ ^(dcf) / dx)(d(j) / dy) . Within the 
Gaussian approximation, the four-body correlation between <fi appearing in the right- 
hand side of Eq. (JT7J) can be decoupled into the product of the pair correlation of (f). 
Noting that the Fourier transform of Eq. (jSJ) gives 0k(^) = 47rZen(e/c 2 )~ 1 '0k(t) and using 
([To]) , we can readily integrate over time in Eq. (fTTj) to find 

k R T , r d 3 k klk 2 



V = Vo + —f^^D 



x "y 



D D J (2tt) 3 k\W + n 2 D y v ' 

Performing the integration over k, we finally obtain 77 = r]o + cikBTno/D, where the 
numerical constant c\ is given by 

1 f 2n • 2 2 , r ■ 5n m f°° x2dx 1 

ci = -, — -r / sin 09 cos ifid(p / sin vdv / ^ = . (19) 

{2nf Jo * * Wo Jo (1 + x 2 ) 3 480tt v ; 

The relative diffusivity D in Eq. (|T3~j) has been originally introduced in the context of 
"collective" meaning. In the limit of the infinite dilution, however, D coincides with 
the self diffusivity of a tagged ion, D s , which is related to the mobility of the tagged 
ion \x via the Einstein relation D s = fikgT. Therefore, by defining the friction constant 
of a tagged ion as 7 = /i -1 , one can recover the limiting law originally derived by 
Falkenhagen for binary electrolytes j3] (and later rederived by Fuoss and Onsager for 
more general multicomponent cases jSllZj), which is written as 

A "» = "-"» = S- < 20 > 

This simple scaling concludes the well-known fact that the viscosity enhancement is 
proportional to the square root of the electrolyte concentration. A similar derivation of 
Eq. (J2UJ) can be also found in the recent literature [121 US] • 
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4. Nonlinear Effects 



Let us proceed to investigate how the electroviscous effect is modified in the presence of 
steady shear flow. In the case of simple shear u = jye x , e x being the unit vector along 
the x axis, the diffusion equation (fT3|) in the Fourier space becomes 



d_ 

dt 







dk 



Mt) = -D{k 2 + 4)^k(*) + W). 



(21) 



Applying the standard manipulation described somewhere 
Eq. (J2T|) in the form 



, we find the solution of 



Ok(s)(t - s)G(k, s)ds, 



where k(t) is the time-dependent wave vector defined by k(i) 
Green's function is given by 

rt 



k + jk x te 



G(k, t) = exp 



-D [ ds(k 
Jo 



(22) 
and the 

(23) 



Using the correlation property Eq. 
given by 



we find that the steady-state structure factor is 



S(k) = 2L / k 2 (s)G 2 (k,s)ds. 
Jo 



(24) 



If we introduce the dimensionless variables q = (D/jY^k and I = js, Eq. (f2*4"|) can be 
rewritten as S(k) = x5*(q), where S(q) is explicitly given by 



S(q) 



1 



dl exp 




— ] I + q x q y l 2 + -qll 3 



(25) 



The Debye relaxation time t d defined by r D = \ 2 D /D = eksT/ (47iZe 2 n t otD) represents 
the time scale at which the deformed ionic atmosphere around an ion recovers its 
equilibrium spherical form. The dimensionless quantity is therefore regarded as 
the control parameter which measures the relative strength of the shear flow (external 
perturbation) to the stability inherent to this system. Numerically evaluated S'(q) for 
two different values of r^j is displayed in Fig. [TJ For td7 > 1, the structure factor 
is considerably deformed by the shear, and rapidly reaches the unity as q x increases, 
reflecting the enhanced homogeneization along the flow direction by the shear. 

These dramatic changes of S'(k) due to shear should also affect the ionic contribution 
to the viscosity. According to Eq. (|12j) . the excess shear stress is given by 



a 



4tt I 



/ AnZen\ 2 



1 



S(k). 



(26) 



e / J (2vr) 3 A; 4 

As is obvious from the spherical symmetry, a xy [^ = 0) is zero. However, to avoid 
apparent divergence of the integrand in Eq. ()26|) in the limit of s = 0, we hereafter 
consider the quantity Acr^^) = o- xy (j) — a xy {f)). Introducing the polar coordinate and 
integrating over k = |k|, we have 



A^(7) = o 



1 ksTn 



8 (2tt) 5 /2 j di/2 



2- 



dip / d9 sin 9 a(9, ip) 
Jo 
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Figure 1. Contour plots in the shear plane of the steady-state structure factor S(q) 
for shear rates (a) td7 — 5.0 and (b) 7737 = 0.5. 



x / dse 
lo 



~Dk 2 d s 



l + 2a(9,(p)rs + b(6,(p)^ 2 s 2 

(s + a(8, p)js 2 + ±6(0, <^)7 2 s 3 ) 3 / 2 s 3 / 2 



.(27) 



where a(9,ip) = sin 2 9 cos <p sin ip and b(9,<p) = sin 2 # cos 2 have been introduced. 
Integrating further over s by parts, we finally arrive at the following result after some 
algebra, 



Ar?(7) 



-Ao- xy (j) 

7 



Ari H(r D j), 



(2f 



where At)q = (480ttD) 1 ksTK,D is the excess viscosity derived in the last section. The 
scaling function H(x) is written as 

h(x) = r 

where 



-2u z 



J(xu 2 )du, 



(29) 



J( y ) = d<p I" d9 sin 9 a(8,<p) 
Jo Jo 



[l + a(9,<p)y + m,ip)y 2 y/i 



(30) 



By expanding J(xu 2 ) in powers of x, one can easily check H(x) = 1 — b\x 2 + • • • for 
i « 1, where b\ = 15/128. The present theory therefore correctly reduces to the 
Falkenhagen-Onsager-Fuoss expression of the excess viscosity in the zero-shear limit. 
For small td7, the excess viscosity is the analytic function of 7 and exhibits a weak 
shear-thinning given by 

I --^r(r D7 ) 2 + •••!• (31) 



A?7(7) = Ar/o 



128 



5. Discussion 



The numerical evaluation of the scaling function H(x) is shown in Fig. As can be 
seen in Fig. |2] (a), the excess viscosity At] rapidly decreases much below the intrinsic 
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Figure 2. (a) Plots of numerically evaluated scaling functions H(x) and xH(x). (b) 
Comparison between (numerically evaluated) H(x) and the fitting function I(x) [Eq. 



value Ar] as the shear rate exceeds 7 C = t d . The overall behavior of H(x) up to x ~ 20 
are fairly well fitted by the empirical function I(x) of the form 

= [rrW' (32) 

where A = 0.17 and a = 0.46 in Fig. 121(b). For small enough x, however, h,2(x) = \-b\x 2 
(the Taylor expansion of H(x) around x = which we have obtained in the last section) 
actually fits better to the numerical date of H(x) than I(x), as shown in the inset of 
Fig. 121(b). Roughly speaking, therefore, the excess viscosity behaves as Ar] ~ 7 _2q for 
1 < T£>7 < 20, and even more strongly decreases as ~ ^~ 5 / 4 for larger Tjj'y. 

Correspondingly, the excess shear stress Aa xy saturates into its maximum value 
Aa m ~ A?7o7 c ~ CiksT/X^ at r D j ~ 10, and then slowly decreases for larger 7. If 
one considers the "effective surface tension" a s of an ionic atmosphere of size Ad on 
analogy with a droplet in an immiscible fluid, one would obtain a s ~ ksT '/X 2 D from the 
energy balance condition. Assuming that the pressure difference between interior and 
exterior of the droplet can be balanced by the imposed shear stress, a s /\o ~ Ao~ m , one 
immediately recovers the above expression of Aa m . 

Let us examine the magnitude of the predicted shear-thinning effect. For the zero- 
shear viscosity at higher concentrations, the following empirical equation: 

At7o(c) = A^d + Bc, (33) 

provides a better description of the experimental results. The constant A was evaluated 
from experimental data for various salts by plotting a suitable form of Eq. (jHHJ), 
and the good agreement with the theory has been obtained for aqueous solutions [3]. 
However, the absolute magnitude of this contribution to the viscosity i^A^fc term) is 
found to be rather small (as can be seen in the large number (4807r) in the denominator 
of Eq. (|20j) ) . Usually, more significant contribution comes from the ion-solvent 
interaction, which is represented by the Be term in Eq. (}3*3*|) . but there is no satisfactory 
theoretical explanation of the coefficient B in Eq. (J33|) so far. The present study 
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suggests that, at strong shear as td7 ~ 10, the coefficient A decreases to 20% of 
its original value. Therefore, while a direct experimental observation of the shear- 
thinning would be hard because of the vanishingly small value of A, the modified 
fitting function At] /r] = 1 + Be, which tells the linear dependence of the viscosity 
upon the concentration, might be best fitted to experimental data under strong shear 
conditions (provided that the linear concentration dependence arising from the ion- 
solvent interaction is not significantly changed by the shear). Although indirect, this 
may serve a possible experimental verification of the shear-thinning effect in future 
experiments. 

At very large shear td7 3> 1, the ionic atmosphere is no longer ellipsoidal because 
it is dramatically stretched by the shear, and dissipates thermally. The considerably 
deformed atmosphere may not support the shear stress any more, and Ao~ xy itself is thus 
decreased. If we use typical values of the diffusion constant D ~ 10 -5 cm 2 sec -1 and the 
Debye screening length A^ ~ 3 x 10 -7 cm for 0.01 M solutions of strong 1:1 electrolytes 
(like NaCl in water), we find 7 C ~ 10 8 sec -1 , which is unrealistically large. This simple 
estimation suggests that, unfortunately, even the strong shear condition td7 ~ 1 might 
be difficult to be realized in ordinary experimental setup. In other words, the ionic 
contribution to the viscosity may not be significantly influenced by the shear in most 
realistic conditions. 

Experimentally, however, a near-critical ionic fluid, in which the charge density 
fluctuation decays slowly, can be prepared by carefully choosing materials and by 
controlling the temperature [221 E31 EH ESj In that case, the correlation decays as 
exp(—Dftk 2 t), where the (renormalized) diffusion constant may be given by Kawasaki- 
Stokes's law Dr ~ kBT/QirrjC, with £ being the correlation length of the critical 
fluctuations [TH1 HHJ - Although a rather crude argument, we naively expect 7 C ~ 3 x 10 4 
sec -1 , which is still very large but is much lower than that for an ordinary electrolyte 
solution, as well as a relatively large magnitude of the electroviscous effect on the order 
of Ar] /r] = ~ 0.1, if we assume i] ~ 40 cP, £ ~ 10 nm, A^ ~ lnm for 

T = 343 K, which are the typical value set for (Bu4NPic)-tridecane-l-ol. (This model 
system shows liquid-liquid phase separation with an upper consolute point near 342K, 
expected to be driven by Coulombic interaction [2H]) Needless to say, in order for a more 
accurate treatment of an ionic fluid near a demixing or a gas-liquid transition, we have 
to take into consideration the density gradient term (Vip) 2 in F which accounts for the 
spatial heterogeneity of critical fluctuations (and becomes more and more important 
as the system approaches the critical point). Our rough estimation therefore might 
be appropriate only in pre-stages of a critical regime. In this connection, the free 
energy Eq.@ including the gradient term is mathematically equivalent to that of 
diblock copolymer systems derived by Ohta and Kawasaki [2B], in which the long-range 
interaction arising from the connectivity of different chemical sequences in a polymer 
chain plays a crucial role. Its close similarity to charged fluids has also been pointed 
out in the literature [2E]. Near-critical properties of that system under shear in one 
phase region was also studied by Onuki based on the same dynamic equations employed 
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here |2Z] • On the other hand, in molecular dynamics simulations, accessible shear rates 
are much larger than those in real experiments, and the strong shear condition will be 
easily realized Then the pronounced shear-thinning behavior predicted here 

would be more likely to be checked by computer simulations. 

The increase and saturation of ionic conductivity of electrolyte solutions at very 
high applied voltage is known as Wien effect EDI EH]- This increase has been 
shown in general to arise from two different effects, namely, a destruction of the ionic 
atmosphere and a modification of the dissociation kinetics 30]. In strong electrolytes 
the former effect can dominate over the latter. It is instructive to remark that the 
similar mechanism works in the system we studied here, that is, the destruction of the 
ionic atmosphere by strong shear renders the increase and saturation of the reciprocal 
of the viscosity, i.e., the fluidity. 

6. Conclusion 

We now summarize our results. Effects of shear flow on the viscosity of electrolyte 
solutions are studied theoretically within the framework of fluctuating hydrodynamics. 
Our calculations confirm the Falkenhagen-Onsager-Fuoss limiting law that the ionic 
effects to the viscosity is proportional to the square root of the concentration of the 
electrolytes in the dilute and the zero-shear limit. We extend this classical result for 
finite shear by deriving the analytic expression of the excess nonlinear shear viscosity. 
The electrostatic contribution to the viscosity is drastically reduced under the strong 
shear condition td7 > 1> which however might be difficult to access experimentally. 
Because the use of Debye-Huckel theory here limits the applicability of our treatment 
to low concentration solutions, the viscoelasticity of concentrated electrolyte solutions 
under shear is an important remaining problem. Studies in this direction using Brownian 
Dynamics simulation is currently in progress. We hope that this work will prompt 
further studies on the nonlinear transport theory of electrolyte solutions in large velocity 
fields, which might constitute the counterpart of the conductivity studies in large driving 
fields. 

Acknowledgments 

The author would like to thank Y. Murayama and M. Sano for valuable discussions. 
This work was supported by Grant-in-Aid for JSPS Fellows for Young Scientists, from 
Ministry of Education, Culture, Sports, Science, and Technology, Japan. 

Appendix A. Derivation of the dynamic equations 

In this Appendix, we describe the derivation of the dynamical equations (JHJ) and (jZJ) 
based on a two-fluid model for electrolyte solutions. We will specifically consider 
strongly dissociative salt solutions such as NaCl and KC1 in water. Such a solution 
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should be modelled as a three-fluid model, namely, as a mixture of dissolved positive 
ions, negative ions and the solvent. For simplicity, however, we will start from the 
two-fluid model for positive and negative ions, while regarding the solvent as the 
structureless background with a static dielectric constant e. A number of approximations 
and simplifications involved at this stage are to be critically checked in future studies. 
(The present formulation is convenient to directly obtain a mutual diffusion equation 
for electrolyte solutions coupled with the Navier-Stokes equation. However, in order to 
obtain a diffusion equation for each ion species with vanishing electrostatic interactions, 
we also have to take into account hydrodynamic drag forces between ions and solvent 
molecules, which are not explicitly included in the present formulation. A more precise 
and systematic treatment for diffusion equations can be found, for example, in Ref. j32j.) 
The number conservation law for the ionic species a = ± is written as 

^ + V-(n n v a ) = 0, (A.l) 
while the total density n tot = n + + n_ obeys the usual continuity equation 

+ V ■ (n tot v) = 0, (A.2) 
where the average velocity v is defined by 

ntotv = + n_v_. (A. 3) 

Assuming that a fluid under consideration shows a vanishingly small compressibility 
in long-time scales, we set V • v = (the incompressibility condition), which implies 
from Eq. (|A.2|) that the total number density is constant, n tot = n tot = (Z + l)n, for 
a sufficiently slow fluid motion where dtn to t can be neglected. Introducing a relative 
velocity between the positive and negative ions by u = v + — v_, we can write the 
velocities of the two ion species as 

v + = v H u, and v = v u. (A. 4) 

n to t ntot 

Since we have defined a scaled charge density fluctuation as Zenifj(r) = Zen + (r) — 
en_(r), we find from Eqs. (|A.1|) and (jA.4|) the equations for ip: 

1 ?iifi_ 



^ + v-(?H = -v 



-u 



(A.5) 



Z n 2 

On the other hand, the equations of motions for the two ions species may be given by 

= -n + Vjz+ + 77+V 2 v+ - Cu, (A.6) 

3~v 

mn_— — = — ra_V/i_ + r]_V 2 v_ + ^u, (A. 7) 

at 

where m is the equal mass for two ion species, rj a is the viscosity, \i a is the appropriately 
defined chemical potential, and ( is the friction constant between the two ion species. 
Note that £ may be given in terms of the friction constant of the each ion species ( a by 
= Cp 1 + (T 1 . Subtracting Eq. (jA.7|) from Eq. (jA.6|) . we obtain the equation for the 
relative velocity u as 

— = V(/i + - - ^ — + — u, (A.8) 

ot m m \n+ n_ 7 
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where the viscosity terms have been neglected. Since we are concerned with very slow 
motion whose characteristic frequencies are much lower than ((l/n + + l/n_), we can 
set d t u = in Eq. (|A.8|) to obtain 

u = -(— + —) iv(/i+-/i_). (A.9) 
\n + n-J Q 

Here the chemical potential difference between two components are connected to the 
free energy as 

SF 5F 5F 
fI+ - fI _ = — __ = C onst + (A.10) 

where F is the free-energy functional in our model, and \ = (Z + l)/(Zn). In 
differentiating F with respect to n + (ri-), we fix n_ {n+) in Eq. (jA.10|) . This relation can 
be directly checked by using the explicit form of F given in the main text. Substituting 
Eq. (1A~9|) into Eq. (lA~5l) with Eq. (jA~Tf)jl . we obtain 

^ + V-(^v) = LV 2 ^(/3F), (A.ll) 

where (3 = l/ksT and the kinetic coefficient 

L = Wxfafi-) 2 = kiF (A12) 

Z C n 2 n tot ( 

is independent of ip. On the other hand, by adding Eq. (jA.6|) and (jA.7|) . the equation 
for the average velocity v is found to become 
d~v 

mn tot — = -(n + VyU + + n_V/i_) + (77+ + r/-)V 2 v, (A.13) 
at 



where the use of Eq. (jA.4|) has been made, and only fluctuations up to the linear order 
have been retained. Again, by using the explicit form of the free energy functional, one 
can directly check the relation 

5F 

n + V/i+ + n_V/i_ = Vp-^V— , (A.14) 

dip 

where p can be regarded as the hydrostatic pressure of the fluid. At a final step, 
adding random thermal noise sources 6{r,t) and f(r, t) which ensure the appropriate 
equilibrium correlations of ip and v to Eqs. (jA.ll|) and (jA.13|) . we obtain a closed set 
of hydrodynamic equations for ip and v 

= _ v . V ^-LV 2 ^(/3F)+fl, (A.15) 

/^7 V = ~ V P ~ k B TijV—((3F) + Vo V\ + f, (A.16) 

at dip 

where p = mn tot is the average mass density, r] = r/ + + is the total viscosity of the 
fluid, and where the incompressibility condition V ■ v = has been used in Eq. (jA.15|) . 
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